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sult shows that the frequently used approximations, for instance to eliminate strongly 
>• coupled optical cavities, are mathematically consistent. 

o 

(N 

O PACS numbers: 42.50.Ct, 03.65.-w, 02.30.Yy, 07.07.Tw 

(N 



> 

X 



a ' Electronic mail: jug@aber.ac.uk 

b ) Electronic mail: 

c ) Electronic mail: 



Hendra.Nurdin 


Sanu.edu.au 


sew08@aber.ac.uk 





1 



I. INTRODUCTION 



Adiabatic elimination is a standard modeling procedure adopted when dealing with sys- 
tems that have both slow and fast variables. Here one considers the limit in which the 
fast variables are effectively relaxed to instantaneous equilibrium values, which may in turn 
depend on external influences, and an effective dynamics may therefore be deduced for the 
slow variables. The problem becomes more involved when the system is driven by stochastic 
influences. In quantum optics, fast oscillators driven by quantum input processes may be 
eliminated from the dynamics using the limiting procedure that they are strongly coupled 
to the input field processes. The first rigorous account of this limit was given by Gough and 
van HandeP^ and the resulting reduced open dynamics for the slow degrees of freedom were 
obtained. Extensions of this result to general nonlinear models with a slow-fast time scale 
separation were given subsequently by Bouten, Silberfarb, and van Handel™. 

Adiabatic approximation is frequently used to simplify the description of a model. In this 
paper we aim to investigate the correctness of applying component- wise adiabatic elimination 
in quantum feedback networks with Markovian components. Here several quantum systems 
may be connected by passing the output noise from one component in as input to another. 
In the zero time delay limit we may model the network as a global Markovian systerrP^. 
For a certain class of quantum networks and under certain conditions, we show that the 
instantaneous feedback limit used to obtain a Markovian quantum feedback network is 
indeed compatible with the component-wise adiabatic elimination procedure. This is the 
ideal situation one would require for modeling purposes, however, the conclusion is not 
immediately obvious when treating individual cases, particularly when the architecture of 
the network becomes complex. We show that for both limits the form of the coefficients 
of the quantum stochastic differential equation (QSDE) describing the limit evolution can 
be formulated as a Schur complement of pre-limit coefficients. Commutativity of the Schur 
complementation procedure then ensures the commutativity of the adiabatic elimination 
and instantaneous feedback limits. 

In section [TT] we shall review the rigorous results that exist for adiabatic elimination of 
oscillator components and adapt the results to deal with multiple oscillator elimination (the 
proof is deferred to the Appendix). We show commutativity of the limits for a simple cascade 
of components and for components in a non-trivial feedback loop. In section [Hi] we present 
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the main features of Schur complementation which we shall need, and show that both limits 
involve Schur complementation procedures. The proof of commutativity of the limits is then 
established in section HVl 



Notation. In this paper we will use the following notation: i denotes y— T, kerX (or 
ker(X)) denotes the kernel of an operator X, imX (or im(X)) denotes the image of an 
operator X. Also, * denotes the operator adjoint. For instance, if X = (Xi,X 2 , . . . ,X m ) is 
a row vector of operators X\, X 2 , . . . , X m on some common Hilbert space then X* is column 
vector given by X* = (X*, X*, X* m ) T . Also, Re c (or Re(c)) and Imc (or Im(c)) denote 
the real and imaginary parts of a complex number c, respectively. 

II. MODELS IN QUANTUM OPTICS 

A. Quantum Input Components 

The standard motivation for quantum stochastic evolutions in physical models has been 
via traveling quantum fields interacting in a Markovian fashion with a given quantum me- 
chanical systenP. The fields may be described by idealized annihilation and creation opera- 
tors bj (t) and bj (t)* respectively, for j = 1, ■ ■ ■ ,n, assumed to satisfy singular commutation 
relations 



These are sometimes referred to as quantum input processes. From these we may define 
regularized operators 



The older, mathematically rigorous approach is that of Hudson and Parthasarathy which 
realizes the open Markov dynamics of a system with Hilbert space f) through a dilation 
to a unitary evolution on a larger space f) £g> Specifically $ is Bose Fock space over 
.8 <g> L 2 [0,oo) where ^ = C n is the colour, or multiplicity, space of the quantum inputs. 
The processes Bj (•) ,Bk (•)* , Aj-& (•) are then realized as concrete creation, annihilation and 
second quantization operators on 

We shall now work in the category of such models: each element of the category will be 
an open quantum system modeling a quantum device. A single component with intrinsic 



[b j (t),b k ( 8 )*] = 6 jk 8(t-s). 
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Hilbert space f) is modeled as an open quantum system with external driving space $ - the 
Bose Fock space over a one-particle space .ft® L 2 (M + ). As above, .ft is the multiplicity space 
of the Bose noise field, and we shall restrict to finite multiplicity for each component (.ft = C n 
for some multiplicity n). Taking {ej}™ =1 to be a fixed orthonormal basis in .ft, we realize 
Bj(t) as the annihilation operator B(ej £g> l[o,t]) on 3, with Bj(t)* the creator. The process 
Ajk(t) is then the differential second quantization of the one-particle operator \ej) (e^l ®l[o,t] 
on .ft® L 2 (IR + ) where l[ 0) t] denotes the operation of multiplication by l[o,t]- We remark that 
we have the continuous tensor product decomposition 

for each t > 0, where #t] is the past noise space (Fock space over .ft £g> L 2 [0, t}) and #[t is the 
future noise space (Fock space over .ft <g> L 2 [t, oo)). A process X(-) on f) £g> # is then said to 
be adapted if for each t > 0, X(t) acts trivially on the future factor #[t- 

The Hudson-Parthasarathy theory of quantum stochastic processes^ gives a non- 
commutative generalization of Ito's stochastic integral calculus. With differentials dBj(t), 
dBl(t), dAjk(t) understood as being Ito incremenlP^ (i. e . ; they are "forward looking", 
dX(t) = X(t + dt) — X(t) where X can be any of Bj, B%, Aj k ), we obtain the following 
quantum Ito table^ for second-order products of the quantum Ito differentials (omitting 
the dependence on t for brevity) 



X 


dB 3 


dAjk 


dB* k dt 


dB l 





SijdBk 


Sikdt 


dAi m 





SmjdAik 


SmkdBf 


dB* m 











dt 












The most general form of a unitary adapted process U(-) on f)®^, with time-independent 
coefficients, will occur as the solution of a quantum stochastic differential equation (QSDE) 
of the form (adopting a summation convention) 

dU(t) = {K®dt- L*S jk <g> dB k (t) + Lj <g> dB^t)* + (S jk - 5 jk ) <g> dA jk (t)} U(t), (1) 

with U(0) = I, and where the damping term is 

K = - l -L)L 3 -iH. (2) 
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We are required to take S = [Sjk] G B(i) <S> A) to be unitary and H to be self-adjoint. The 
operators Lj and H are assumed to have a common dense domain in f), which holds in 
particular when these operators are bounded. In the case that they are unbounded, they 



will be of a particular form which will be given in (10). 



output 



-o (S,L,K) 



input 



FIG. 1. Single component 

From our point of view the category of all possible components is parameterized by t),n 
and the possible triples (S,L,K) as above. It is convenient to collect all the coefficients in 
the QSDE Q into a single operator G G B(Sj) where 

£ = &®(C©j!). (3) 

With respect to the decomposition C © ^ we specifically define G to be 

K -L*S 



G 



(4) 



L S-I 

In this representation, G appears as a (1 + n)-dimensional square matrix with entries in 
5(f)). 

In Fig. [T] we sketch the open system as an input-output device specified by the triple of 
operators (S, L, K). The output fields is defined to be the canonical processes 



B°"\t) = U(ty[I®B j (t)p(t)- 



(5) 



B. Systems in Series 

Let us consider a pair of systems (Sj,Lj,Kj), j = 1,2, connected in series as shown 
in Fig. [2j (Note that we technically do not require the observables of the two systems to 
commute!). 
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-o(S 2 ,L y igo 



- CS^L^KJo- 



FIG. 2. Systems in series 

In the instantaneous feedforward limit, the pair can be viewed as the single system shown 
in Fig. wlwith overall parameters^ 



(S scv ,L sev ,K ser ) = (S 2 ,L 2 ,K 2 ) < (S l ,L 1 ,K 1 ) (6) 

where the series product < is the associative (though non-commutative) product given by 
the explicit identification 

Sser = S2S1, (7) 
-^ser = Li2 + S%L\, (8) 

K„ 3t = K 1 +K 2 -L* 2 S 2 L l . (9) 



O 2 o- 



input 



I.F. 



output 



Series 1 & 2 



3- 

— ' input 



FIG. 3. Systems in series: the upper setup describes two systems connected in series with a time 
lag r > in the interconnection from system 1 to 2. In the instantaneous feedforward (I.F.) limit 
we consider r — > in which case we obtain an effective single component model again. 

Note that if Hj (j = 1,2) are the Hamiltonians of the separate systems then the damping 
operators are Kj = —^L*Lj — iHj and the effective Hamiltonian in series is then given by 

H aeI = H 1 + tf a + Im(LS£ 2 £i). 
C. Adiabatic Elimination of Oscillators 



We suppose that the system consists of local oscillators having Hilbert space f) osc and 
that the remaining degrees of freedom live on an auxiliary space f). The overall Hilbert 



space of the system is then f) ® () osc and we consider an open model described by the triple 
of operators 

S(k) =S®I, 

L(k) = k s ^C j ®a j + G®I, 

j 

K(k) = k 2 ^2 A Ji ® a *j a i + k Zj ®a* + k Xj ®a,j + R®I, (10) 
ji j j 

where k is a positive scaling parameter and S,Cj,G, Aji, Xj, Zj, R are bounded operators 
on f) with A = [Aji] boundedly invertible. Here cij is the annihilator corresponding to the 
j-th local oscillator, say with j — 1, • • • ,m. 

As /c — )• oo the oscillators become increasingly strongly coupled to the driving noise field 
and in this limit we would like to consider them as being permanently relaxed to their joint 
ground state. The oscillators are then the fast degrees of freedom of the system, with the 
auxiliary space f) describing the slow degrees. In the adiabatic elimination fc->oowe desire 
a reduced description of an open system involving the operators of f) only, with the fast 
oscillators being eliminated, as illustrated in Fig. |4} The ground state for the ensemble of 
m oscillators will be denoted by |0) osc . 

Define X to be the row vector of operators X = X 2 , . . . , X m ) and Z to be the column 
vector of operators Z = {Z\, Zi, ■ ■ ■ , Z m ) T . Also, we say that a matrix M = [M ? 7] :7 y =lr .. jm , 
with Mji bounded operators on f), is strictly Hurwitz stable if Re(0, M<j>) < for all ^ (p G 
f) eg) C m . Then we say that an open Markov quantum system with parameters of the form 



(10) is strictly Hurwitz stable if the matrix [Aji] is strictly Hurwitz stable. We first have the 



following result: 



Theorem 1 Let U(-, k) be the unitary adapted evolution associated with the triple (S(k), L(k), K(k)) 



appearing in (10), and define the slow space as l) s = f) <8> {C |0) osc } . If the operator 
Y = Ylji Aji <8> a*ai has kernel space equal to the slow space, then we have the limit 



lim sup 

fc^oo o<t<T 



U(t,k)$ - U(t)$ 



0. 



for all T > and $ G f) s <S> 3, where U(-) is the unitary evolution associated with the triple 
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S® |0>(0| OSC ,L® \0)(0\ OSC ,K® |0)(0| osc and 



S=(I + CA- l C*)S, 
L = G — CA- l Z, 
K = R- XA~ l Z. 



(11) 



Remark 2 For ease of notation, we will drop the factor "• <g> |0)(0| OS c" it is obvious that 
in the limit we are always relaxed into the fast oscillator ground states. Therefore we can 
simply think of the limit QSDE as having initial space f) and coefficients (S,L,K). 



Remark 3 A sufficient, though not necessary, condition for the kernel space of Y to equal 



the slow space is that the matrix A be strictly Hurwitz, see Lemma 15 



The result is a generalization of what has been established for the single mode case^ 
where the main result is stated for weak convergence of the unitaries, but this automatically 
extends to the strong convergence above. There the techniques were based on a quantum 
central limit theorempl which have been shown to extend to the multimode situational. We 
shall give a proof the theorem in the Appendix, exploiting the theory of singular perturbation 
of QSDEs developed by Bouten, van Handel, and SilberfartP^. 

In the following, we shall drop the tensor product symbol for notational convenience. 
Furthermore we shall introduce the vectorial multi-mode notation 



We therefore write simply 



S(k) = S, L(k) = kCa + G, K(k) = k 2 a*Aa + ka*Z + kXa + R. 



If we take the Hamiltonian to be 



H(k) = k 2 a*na + ka*Y + kY*a + 6 



Auxiliary System 
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A.E. 



oscillator 



output 



o Auxiliary System o 



input 



output 



input 



FIG. 4. The setup on the left shows a system of an oscillator and auxiliary component with 
the oscillator coupled to a quantum field input. In the limit where the coupling of the oscillator 
becomes infinitely strong, we may adiabatically eliminate the oscillator to obtain an input acting 
directly on the auxiliary component. This is sketched in the setup on the right. 



then 



A 
Z 
X 
R 



-c*c - in, 



--c*G-ir, 

— G*G - iQ. 
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In particular we note the identities 



A + A* = -C*C, 
X + Z* = -G*C, 
R + R* = -G*G. 



(12) 



We present a naive derivation of the limit form appearing in Theorem [TJ with the proof 
presented in the Appendix. In the interaction picture we have the quantum Langevin equa- 
tion 

a = l -L(k)* [a, L{k)\ + ^L(k)* [a, L(k)) - i [a, H(k)) - [L(k)*S, a] b m 



1 



1 



= ~k\-C*C + iti)a + k(--C*G - iT) - hC*Sb- m , 

where b m is an input quantum process satisfying [b in (t), b in (s)*] = S(t — s). Likewise the 
input-output relations are 



-'out 



Sb in + L(k) = Sb in + (kCa + G), 



where b out is the output quantum white noise field. 

We note that we may rewrite the Langevin equation as -ra = —kAa + Z — C*Sb- m and 
one argues that as k — > oo the left hand side vanishes, so that the right hand side may be 
rearranged as 

ka^ A-\C*Sb in - Z). 

The common interpretation of this is that the (scaled) oscillator mode becomes "slaved" 
to the input field: usually this argument is given with k fixed to unity and while clearly 
mathematically problematic nevertheless, rather miraculously, yields the correct answer. 
Making this substitution in the output relations, we reasonably expect that 



Km = (I + CA- x C*)Sb- m + {G — CA~ l Z) 
= Sb in + L. 



This justifies the form of S and L. The form of K may be deduced by substituting ka ~ 
A~ 1 (C*Sb- in — Z), and the conjugate relation, into the Langevin equation for any operator 
acting nontrivially only on the space f). (There is a potential operator ordering ambiguity 
here, and the appropriate choice is to substitute ka and ka* in Wick ordered form!) 



D. Adiabatic Elimination and Systems in Series 

The aim of this section is to determine whether the limits of adiabatic elimination and 
instantaneous feedforward do in fact commute, as illustrated in Fig. |5j While this is often 
assumed in quantum optics models, it is certainly far from obvious. At this stage, however, 
we are able to reduce the question to a direct computation. 

Let us represent the local oscillators a\ and a2 in a combined manner as 



a i 
a 2 
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Component 2 Component 1 /^E. 



output 



2 o 



© 1 



input 



output 



input 



I.F. 



output 



Component 1 & 2 


< 





A.E. 



V output >— 



input 



input 



FIG. 5. The picture illustrates the main result that we which to prove here, namely that the 
adiabatic elimination (A.E.) and the instantaneous feedforward (I.F.) limits can be interchanged. 

Then the first system is to be represented as (S 1 (k), Li(k), Ki(k)) where 
Si(k) = S 1 , 

Li(ife) = fcCiai + G ± = k[C u 0]a + G u 
K^k) = k 2 a\A x a x + ka* l Z 1 + kXa x + R x 

+ k [Xi, 0] a + R 1 . 





A x 




Zi 


k a 


a + ka* 













Likewise, the second system is then represented as 
S 2 (k) = S 2 , 



L 2 (k) 


= k[0,C 2 ]a 


+ G 


2, 



















K 2 (k) 


= k 2 a* 






a + ka* 









A 2 




z 2 



+ k[0,X 2 ]a + R 2 . 



1. Adiabatic Elimination followed by Instantaneous Feedforward 

If we perform the adiabatic elimination first then we arrive at the two systems (j 

Sj^il + CjAj^Sj 
L i = Gj - CjAj 1 Z j , 



1,2) 



1 ' / _ l Zj- 



The instantaneous feedforward limit is then given by the series product 
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S — S 2 Sxj 
L = L 2 + S 2 Lx 
K = k x + K 2 -L\S 2 L 1 . 



(13) 



2. Instantaneous Feedfoward followed by Adiabatic Elimination 

We perform the series product (S 2 (k), L 2 (k), K 2 {k)) < (Sx(k) , Lx(k) , Kx(k)) first to ob- 
tain (S 8ei (k) , L scr (k) , K scv (k)) where 

S S er(k) = S 2 Si, 

L BeI (k) = L 2 (k) + S 2 L!(k) 

= k[3 2 C 1 ,C 2 ]a + G 1 + S 2 G 1 , 
K SCI (k) = K x {k) + K 2 (k) - L^kyS^L^k) 
A x 
-C 2 S 2 Cx A 2 



2 * 



k z a 



a + ka* 



Zx 

Z 2 — C 2 S 2 Gx 



+k [Xx — G 2 S 2 Cx, X 2 ] a + Rx + R 2 — G* 2 S 2 Gx- 

Now, adiabatically eliminating the oscillators leads to the effective model (S seI , L ser , H se r)- 
Here we have 



S ser ={I+[S 2 Cx,C 2 ] 



Ax 
-C 2 S 2 Cx A 2 



-1 r 



CIS* 



c 2 



)S 2 S_ 



1- 



L S er = (Gx + S 2 Gx) — [S 2 Cx, C 2 \ 

Kser = (Rl + R2 — G* 2 S 2 Gx) 

— [Xx — G* 2 S 2 Cx-,X 2 ] 



Ax 
-C 2 S 2 Cx A 2 



Ax 
C 2 S 2 Cx A 2 



Zx 

Z 2 — C 2 S 2 Gx 



Zx 

Z 2 — C 2 S 2 Gx 



(14) 



E. Commutativity of the Limits: Systems in Series 



The matrix inverse appearing in ( 14 ) is easily computed as a special case of the well-known 



formula for the inverse of block matrices (the earliest reference is credited to Banachiewic^, 
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see subsection HI A however, like many matrix identities the origins may be considerably 
older) 



A x 
— C 2 S 2 Cx A2 

This yields the explicit form 



-1 -1 



A^ 1 
A^C^CxAj 1 A 2 l 



Sser — (I + S%CxAx 1 C*S 2 + C2A2 1 C 2 S 2 CxA 1 l CiS 2 + C 2 A 2 1 C 2 )S 2 Sx 
— (I + C 2 A2 1 C* 2 )S 2 (I + CxA^C*)Sx 
= S 2 Sx- 



The coupling operator is 



L se r — (Gx + S 2 Gx) — S 2 CxA 1 1 Z 1 

—C 2 A 2 l C 2 S 2 CxA l 1 Zx — C 2 A 2 1 Z 2 + C 2 A 2 l C 2 S 2 Gx 
= (G 2 - C 2 A~ 1 Z 2 ) + (/ + C 2 A 2 1 C*)S 2 (G' 1 - CxA^Zx) 
= L 2 + S 2 Lx- 

Finally we see that 

K ser = Rx + R 2 — G 2 S 2 Gx — X 1 A 1 1 Z 1 

J rG 2 S 2 CxA 1 1 Zx — X 2 A 2 1 Z 2 + X 2 A 2 1 C 2 S 2 (Gx — CxZx, 

We would like to show that this equals Kx + K 2 — L 2 S 2 Lx, now we have 



(15) 



Kx + K 2 - L* 2 S 2 Lx =Rx~ XxA^Zx + R 2 - X 2 A 2 l Z 2 



- (G* - z;a 2 l *c* 2 )(i + C 2 A 2 'C;)S 2 (Gx - CxA^ l Z 



and to compute this we note that A 2 = —\C 2 C 2 — iQ 2 so that 

A 2 U C* 2 (I + C 2 A 2 l C*) = A 2 U (I + ClC 2 A 2 v )Cl 

= a 2 u (a 2 + c;c 2 )a 2 1 c; 

— A~ 1 *(— A*\ A~ x r* 

— — A~ x r* 

— ^2 °2i 



(16) 
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this yields 



K 1 + K 2 - L\S 2 U = Ri- X 1 A l 1 Z 1 + R 2 - X 2 A 2 1 Z 2 

-G* 2 (I + C 2 A 2 1 C* 2 )S 2 {G l - C 1 A^ 1 Z 1 ] 
— Z 2 A 2 ^ C 2 S 2 (G\ — C\A^ ^Z\j. 



We therefore find that 



K scv - (#! + K 2 - L 2 S 2 Li) = {X 2 + G* 2 C 2 + Z* 2 } A?C%S 2 {G X - C X Z 



however this vanishes identically by the second of identities (12) 



We therefore conclude that the model parameters in ( 13 ) are identical with those in ( 14 ) 



therefore the adiabatic elimination and instantaneous feedforward limit commute. 



F. Adiabatic Elimination: In-Loop Device 

Next we want to extend our investigation to situations where we have a non-trivial feed- 
back arrangement as illustrated in Fig. [6] The question again is whether the two limits 
commute. 



Output 



System 

(So, L , H ) 



Input 




FIG. 6. General feedback arrangement: The four port device interacts with one external input, 
producing one external output and interacts with an in-loop device by one internal in- and output 
field respectively. 



We start off with a simple model in-loop, taking the 4-port device to be a beam splitter, 
modeled by a unitary matrix T = [Tj{\ with complex entries and where coupling operators 
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and the systems Hamiltonian are zero, L\ = L2 = Hi = 0. We parameterize the in-loop 
device as 



S (k) = S 
L {k) = kyfya 
K (k) = -^k 2 a* -fa . 



(17) 



and fix the beam splitter (scattering) matrix as (with a real) 



T 



a 



(18) 



-a: 



Output 

Beam Splitter 




Oscillator 



FIG. 7. Oscillator in-loop 

Thus, the in-loop system consists only of a single oscillator coupled to the in-loop 
field and no additional modes coupled to the oscillator. In terms of operator parameters 



Sq, C , G , A , Z , X q , R , see equation (10) 



S (k) = S ®I 

L (k) = kCo®a + Go® I 

K (k) = k 2 A ® a* a + kZ <g> a* + kX ® a + R ® / 
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acting on the space f) sys <g> f) osc , we see that 



So — So 
Ao = ~\l 

c = Vi 

Zq — Xq = Rq 



G n = 0. 



The coefficients for the single input single output device after taking the instantaneous 
feedback limit of the arrangement of Fig. [6] were derived by Gough and James^ and are 
given by 

Sred = Tii + T 12 Sq(I — T22S0) 1 T 21 

-^rcd = T 12 (I — T22S0) L 

H ied = K Q - L* S (I - T 2 2S )- 1 L . (19) 



For the model (17), the limit coefficients after taking the adiabatic elimination limit (see 
Theorem [T| are given by 

So = (l + CqA q 1 Cq) So = —So, 
Lo — Go — CqA q Z = 0, 

K = Ro ~ XoA^Zo = 0. (20) 



Substituting into ( 19 ) we find that the reduced coefficients after the instantaneous feedback 



limit for the model (20) are 



S = a + \/l — a 2 (—So) 



-a)(-S Q 



a — So 
1 — aSo 



(21) 



L = 
K = 0. 



We now exchange the order in which we perform the limits. The instantaneous feedback 
limit of the model before taking the adiabatic elimination limit yields: 

S(k) = Tn + T 12 S (I - T22S0Y 1 T21 = a + (1 - a 2 ) S —^—; 

1 + abo 



L{k) = kVl-a 2 



1 



1 + aS 



K(k) = K (k) - L (ky So ^ 22 L Q (k) = k 2 al [-U + 7t^V 

1 — 00-/22 V 2 1 + ab 



a . 
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where the operator parameters are 



= 7 1 - aS 
2 1 • aS 

= yi - 

1 + aS'o 
, (l-« 2 )^ 
1 + «o 
G = X = Z = i? = 



The A.E. of the I.F. limit model then results in (here |1 + aSo| 2 = (1 + aSg)(l + aSo)) 
S=(l- (!-« 2 )7 2 1 + ^0 ^ / , (1 - a 2 )S c 



1 + aS'op 7 1 — aSo J \ 1 + aSo 
aS£ - 1) (1 + aS ) ( a + S 



(1 + aS£) (1 - aS ) \l + aS c 
a — So 



a Sq 



(22) 



We see that the limits do in fact commute since we obtain the same operator S in (21) 



and (22), likewise for the operators L and K. The apparently miraculous agreement comes 
as a general feature that will be observed in more complex networks. Our approach will be 
to encode both these limits as instances of a Schur complement operation: the miraculous 
agreements that one encounters in a case-by-case study are in fact just by-product of these 
operations. 

If So = 1 then in quantum optics the system (S (k), L (k), K (k)) represents a one-sided 
optical cavity in which the coupling coefficient ky/j of the partially transmitting cavity 
mirror is large (for large k). The calculations of this section show that for large k the 
network in Fig. [7] can be consistently approximated by an effective device that shifts the 
phase of the output field with respect to the input field by an amount determined by the 
parameters of the cavity and beam splitter. Alternatively, one can also think of the network 
of Fig. [7] as approximately implementing a phase shifting device. 



III. ADIABATIC ELIMINATION WITHIN GENERAL NETWORKS 

The situation of two systems in cascade, as depicted in Fig. [2j is the simplest form of a 
nontrivial quantum feedback network. We remark that at no stage of the calculations did 
we assume that the operators describing the first system commuted with those of the second 
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system. Indeed, the series product is valid even if we do not assume that we are dealing 
with separate cascaded systems and is applicable to the problem of feedback into the same 
system. 

In Fig. [8] we describe a somewhat more engorged quantum feedback network featuring 
feedback and feedforward interconnections. For each component of the network, we will 
have the same multiplicity for the input fields as the output fields, though we split up the 
inputs and outputs geometrically to indicate different physical connections for the fields. The 
unitary S for a given component now additionally implies that we can use the component 
to mix the input fields, with a beam-splitter being the very special case where the entries of 
S are just complex constants. This feature introduces the possibility of topologically non- 
trivial feedback loops that were not present in the simple situations of direct feedforward or 
feedback occurring for systems in series. 

We now aim to investigate the procedure of adiabatic elimination of fast degrees of free- 
dom from components in a general quantum feedback network and, in particular, answer the 
question of whether this commutes with the Markov limit in which we take vanishing time 
lags for the various internal fields in the network. The adiabatic elimination of oscillators 
for components in series will then be a very specific case of this general theory. 




FIG. 8. Quantum feedback network 



The essentially mathematical element in the investigation will be that both the adia- 
batic elimination limit and the instantaneous feedback limit for a general quantum feedback 
network are actually instances of a Schur complement of the Ito matrix G. 



A. The Schur Complement 



We begin by recalling some of the basic definitions and notations relevant for Schur 
complements. For general reviews, see the survey article by Oullette^ or the book chapter 
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by Horn and ZhangP. We shall elaborate on several of the well-known results presented in 
the reviews, largely to take account of the fact that we are dealing with block-partitioned 
matrices with operator entries. In particular we give some minor technical extensions where 
we are explicit about the domains, kernel spaces and image spaces on which the operators 
act. 

A B 



C D 



relative to the mxm 



The Schur complement of an (n + m) square matrix M - 

sub-block A is defined to be 

M/A — D — CA- X B 

under the assumption that A is invertible. We note the following elementary formula, due 
to BanachiewicJ^, for invertible M 



A' 1 + A- l B (M/A)' 1 CA- 1 -A~ l B (M/A)' 1 
-{M/A)' 1 CA- 1 (M/A)' 1 



Definition 4 A matrix M is a generalized inverse for a square matrix M if we have 

A B 



A 


B 


-i 


C 


D 





MM M = M. The generalized Schur complement of M 

M/A — D — CA~B. 



C D 



is then defined to be 
(23) 



Lemma 5 The generalized Schur complement M/A is well-defined and independent of the 
choice of the generalized inverse A~ so long as we have the following inclusions of image 
spaces imB C imA and kernel spaces ker A C kerC. 

Note that imB C imA occurs if and only if ker A* C ker B*. (Recall that the image, or 
column space, of a matrix is the span of its columns, or more generally im (M) = ker (M*) .) 

Lemma 6 For two matrices M and N and some generalized inverse M~ of M we have that 

MM .V .V // im .V c im M 

and 

NM~M = N if ker M C ker N. 
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For the proofs of these lemmata, see Horn and Zhangj^ they are a straightforward conse- 
quence of the definition of a generalized inverse and the postulated image/kernel inclusions. 

The Schur complement and Lemmata [5] and [6] above may be generalized to matrices with 
operator entries. Let M be a bounded invertible operator on a Hilbert space and let us fix 
a decomposition fj = for some finite index set 3. We denote by Xj the component 

of a vector x G in Sjj, and Mji the block component of M mapping from Sjj to fj/. For 
A = {a±, • • • , a n } and B = {b±, • • • , 6 m } non-empty subsets of 3 we write 





x a\ 




M aibl ■ 


■ M a , b 


X A = 




,M A}B = 












M a b ■ 


■ M a b 



The single equation Mx = u then corresponds to the coarsest block form My^x? = uj. In 
contrast, the full system of equations J2ie3 ^ji x i = u j gi yes t ne finest block form. More 
generally, we may examine intermediate partitions. Let A and B be non-trivial (i.e., non- 
empty, proper) subsets of 3 then the equation Mx = u may be written as 



M AB M AB , 




x B 




U A 


M A , )B M A ,, BI 




X B i 




U A > 



(24) 



where A' denotes the complement of set A in 3, and the inverse relation is 



X B 




(M- 1 )^ (M~ 1 )b,a' 




ua 


X B > 




{M~ l ) B , A {M~ l ) B i >A > 




u A > 



We now recall the definition of the generalized Schur complement, sometimes also known as 
the shorted operator, in the case where M need not be invertible. 

Definition 7 Let A and B be non-trivial subsets of the index 3, and let C be a non-trivial 
subset of A, and D be a nontrivial subset of B. Furthermore take \A\ = \B\ and \C\ = \D\. 
Suppose that the sub-block M c ,d possesses a generalized inverse denoted by (M c ,d)~ , then 
the Schur complement of M A>B relative to M c ,d is defined to be 

M ab /M c ,d = M A/CtB/D - M a/ cA m c,d)~M C)B/d . 

In the special case where A = B = Z, we shall simply write M/M c ,d for M^^/M^d- 
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The generalized Schur complement is well-defined and independent of the choice of gen- 
eralized inverse so long as the column space ud.(Mc > b/d) is contained in im.(Mc jD ), and 
ker(Mc£>) is contained in ke^M^/c^). In particular, if the conditions of the Lemma [5] are 
met then we may fix a particular generalized inverse such as the Moore- Penrose inverse. 
We also remark that we may readily extend the above notation to the case where different 
direct-sum decompositions of fj are used for the columns and rows of M. 

We shall also require the extension of the Banachiewicz formula to generalized inverses. 
The proof for finite rank matrix operators is due to Marsaglia and StyarP, and may be found 
as Theorem 4.6 in the article of Ouellette^. In the next lemma, we strengthen this to deal 
with general Hilbert space operators. 



Lemma 8 Let M be partitioned according to M 



A B 
C D 



We suppose that im B C im A, 



kerA C kerC, and therefore the generalized Schur complement X = M/A = D — CA~B 
is well-defined and independent of the choice of generalized inverse A~ to A. Then the 
generalized inverse of M is given by 



M~ 



A- + A~BX-CA- -A~BX- 
-X-CA- X- 



Proof. We multiply out MM M in block form. The top left block will be 

(MM-M) n = AA-A + AA~BX~CA-A - BX~CA~A 
-AA'BX'C + BX~C 
= AA~ A + AA~BX~C(A~A - 1) - BX-C^A'A - 1) 
= AA-A + (AA~ - l)BX-C(A7A - 1) 
= A, 

where the last step follows because AA~ A = A and (AA~ — 1)B = under the assumptions 
that imB C imA Similarly 

(MM-M) 12 = B + (AA~ - 1)B + 

(AA~ - l)BX~ C A~ B(AA- - l)BX~D 
= B, 
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since under the assumption imB C imA we have that A A B = B and so (AA —1)5 = 0; 

{MM~M) 21 = C + C(A~A - 1) + (CA~B - D)X~C{A'A - 1) 
= C, 

because of the assumption kerA C kerC we have CA~A = C and C(A~A — 1) = 0; and 

{MM~M) 22 = D-(D- CA~B) — (D — CA~B)X~(D - CA~B) 
= D-X + XX- X 
= D, 

since X = M/A = D — CA~B and XX~X = X. Collecting these results we have that 
MM~M = M, as required. ■ 

Now, as a corollary to Lemma [8] we obtain the generalized Banachiewicz formula: 

m b,a = m a,b + M^ B M A>B ,(M/M AiB )-M AI>B Ml B , 
M BA , = -MX tB M A ^(M/M AB )-, 
M~, A = -(M/M A>B )-M AljB M A>B , 
M B , A , = (M/M AjB )~. 

We now wish to establish the property of commutativity of successive Schur complemen- 
tation as this shall be the main technical result required in this paper. 

Lemma 9 (Successive complementation rule) Suppose that A,B,C is a partition of 
the index set 3 (that is, A,B,C are disjoint non-empty subsets whose union is Z) then, 
whenever the generalized Schur complements are well-defined, we have the rule 

M/M Bu c, B uc = {M/M c ,c)/{M/M c ,c)b,b 

= {M/M B , B )/((M/M B>B ))c,c- (26) 



For the case of matrices over a field where the inverses exist, the first equality in (26) is 
an instance of the Crabtree-Haynsworth quotient formula 6 . The extension of the quotient 
formula to generalized inverses for matrices over a field was given by Carson, Haynsworth 
and MarkhamP^ under some rank conditions, see Theorem 4.8 in the review by Ouellete 11 . 
However, since here we are dealing with matrices with Hilbert space operator entries rather 
than just matrices over a field, we need to extend this result accordingly. To this end, below 
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we independently prove a generalization of the algebraic content of the theorem to matrices 
with Hilbert space operator entries, modulo the conditions for these Schur complements to 



be well-defined which we defer to Lemma 17 in the Appendix. 



Proof. Assume that the conditions of Lemma [17] are in place. Let us first compute 

(M/M CiC )/(M/M c ,c)b,b- 



M/M c ,c = 



M A;A M A>B M A>C 

M BA M B)B M B , C 

M c ,a M C , B M c ,c 

M A;A M \.n 

M BA M BB 



/M c ,c 









M BtC 



(M c ,c)- [M C , A M c ,b] 



= [M L , R - M L , c (Mc,crM c , R } RL&{A B} 
so a second Schur complementation leads to 

(M/M Ci c)/(M/M c ,c)b,b = M AA - M A c(Mc,c)~M c , A 

- (M AjB - M A , c {M CjC )-M C)B )»{M B)A - M B>c {M CjC )-M CA ), 

where we write 5 = (M BjB — M B) c(Mc,c)~ Mc\ B )~ for shorthand. We then compute 

M/M BuC , Bu c 



M/M_ 



BUC,BUC 



M AA M AB M A>C 
M BA M BB M B , C 
M c ,a M C<B M CtC 



M AA - 



M B>B M Bj c 
M c ,b M CyC 

U.x.r, M A 



c 



M B)B M B , C 
M c ,b M Cfi 





M BA 




M C ,A 



however, the inverse can be computed explicitly using the Banachiewicz formula as 

S -EM B ,c(M c ,c)~ 
-(Mc, c yMc, B E (M c ,c)- + (M CiC )-M CiB SM BiC (M CjC )- 

Multiplying out the block matrix readily leads to the same expression already obtained 



for (M/Mc,c)/{M/Mc,c)b,b- The second equality in (26) follows from Lemma 17 and by 
interchanging B and C. ■ 
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B. Instantaneous Feedback Limit as Schur Complement 





Si 





... o 




u 







s 2 


... o 






s = 






'•. 














o s c 







Suppose that we are given a collection of components which have separate descriptions 
(Sj, Lj, Kj) for j — 1, 2, ■ • • , c. We may collect them into a single model (S, L, K) given by 



This just describes the open-loop system where no connections have been made and all input 
and output fields are therefore external. 

To obtain the closed loop description we have to give a list of which outputs are to be fed 
back in as inputs. Algebraically, this comes down to assembling a total multiplicity space 
& = ®j =1 &j and a joint system space f) = <8)S =1 f)j. In this way we obtain a network matrix 
G on f) © (C © &) from the component matrices on [)j <g> (C © 

Once the connections have been specified, we arrive at a decomposition 

where ^ e lists all the external fields and ^ lists all the internal fields. This decomposition 
induces a decomposition of f) as 

[) © (C © £) = [fj © (C © £ c ) ] © ft © #i] 

and with respect to this decomposition, the Ito matrix may be partitioned based on internal 
('i') and external ('e') components as (see Gough and James^ for details) 



G 



^"ee ^ei 



(27) 



In Fig. [9] we sketch the picture that emerges when we subsume all the external fields together 
and all the internal fields together as single channels. 

In the instantaneous feedback limit we find that the reduced model is described by the 
Ito matrix G B(i) © (C © 8^)) given by the Schur complement 



J-"G — G cc — G C i(Gji) 1 G ? ic, 
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(28) 



external 
output 

< 



external 
input 



internal output-input 
FIG. 9. Feedback situation 



related model matrix U 



rather than G and the corresponding feedback 



provided that G\\ exists. We remark that the original version of this formula^ involved the 

K -L*S 
L S 

reduction map was the fractional linear transformation J-U = U ee + U ei (l — Ua)^ 1 !!^. In 
both cases, the condition that G\\ = Sa — I be invertible is necessary for the feedback network 
to be well-posed. 

Remark 10 We note that models studied here all satisfy a Hurwitz stability condition, 
though not necessarily in the strict sense. In general, the feedback reduction need not pre- 
serve the strictly Hurwitz property, and we may obtain conditionally stable modes through 
interconnection. 



C. Adiabatic Elimination as Schur Complement 

We now re-examine the adiabatic elimination of oscillators. For finite k we consider the 
ltd matrix 

K (k) -L(k)*S 
L{k) S-I 



G(k) 



where we write the scaled operators (10) as 



K (k) = [I, ka*} 
L (k) = [I, ka*} 
S = [I, ka*} 



R 


X 




I 


Z 


A 




ka 


G 


C 




I 










ka 



S 




I 







ka 
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Recalling Remark |2| , it is now convenient to use the decomposition f) = f) © f)/ (here f) / 
denotes the subspace of the fast oscillator modes) to write 

\) © (C © £) = \h © (C © J*)l © [f) f © (C © £)] 

and with respect to this decomposition we may now write 

G (k) = k 2 a*gfja + ka*gf s + kg sf a + g ss 



[I, ka*} 



and we set 



9ss g s f 

9fs 9ff 

9ss 9sf 

9fs 9ff 



I 

ka 



(29) 



(30) 



It is easy to see that g is given by 





R 


-G*S 




X 





9ss = 


G 


S-I 


,9sf = 


c 







Z 


-c*s 




A 





9fs = 








,9ff = 









The Ito matrix corresponding to the limit operators (S,L,K) in (11) is then 





k 


-L*S 




R- 


XA- L Z 


-G*S + XA- 1 C*S 


G = 
















L 


s 




G- 


CA- l Z 


S + CA- 1 C*S 



where we use the identity —L*S = — G*S + XA C*S in the upper right corner which relies 



on the trick (16) along with the identities (12). We then observe that 

G = g ss - g sf (g ff )~ g fs = g/g ff 
which is the generalized Schur complement based on the Moore-Penrose inverse 



A 







A' 1 














Indeed, given the specific form here we see from the remarks after the Definition [7] that 
any generalized inverse may be used here. We may then define the adiabatic elimination 
operator as A : G (k) i— )• G = g/gff- 
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IV. COMMUTATIVITY OF THE LIMITS IN GENERAL NETWORKS 



Our first step is to see how the instantaneous feedback limit sits with the adiabatic limit 
starting from a general model with fast oscillators and internal connections which we wish 
to eliminate. 



We have seen from (29 ) that the Ito matrix G (k) may be written as G (k) = [I, ka*]g 



ka 



with g given by (30). suppose that the input fields can be partitioned into internal and 



external fields that corresponds to a partitioning of S as 

^ee Sei 



s 



where Sa is a square matrix pertaining to the scattering of the internal fields to themselves, 
S ee is a square matrix pertaining to the scattering of the external fields to themselves, 
while S e i and Si e pertains to a scattering of internal fields to external fields, and vice-versa, 
respectively. We also partition C and G accordingly as 









G c 


c = 


;G = 




a 







If we wish to decompose this with respect to the external and internal field labels, then we 
obtain 



G(*) = 



Gee (k) G ei (k) 
G ie (k) G n (k) 



similar to (27). We note that these blocks will necessarily have the following structure 



Gee (k) 


= [I, ka*]g 


G ci (k) 


= [/, ka*]g 


G ic (k) 


= 9ee 


I 

ka 


Gn (k) 


= 9n, 





I 

ka 
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with 





Ri 











9ee = 


G l 


s^- 1 


; 9ei = 


a 







z x 


-C*Si 




A 





9ie = 






; ga = 























and 



S e = 


'S'ee 








R — G*S e 




-G*Si 








, Ri — 






'S'ie 








G c S cc / 







G x 



Gi S- u 



Z -C*S 



X 

c e 



I 

ka 



We therefore obtain the feedback reduction 

TG(k) = G cc (k) = [I,ka*} (g/g. 

Conversely, the adiabatic elimination corresponds to 

AG (A;) = g/g, f . 

In this way we see that the essential action is a Schur complementation of the object g either 
with respect to labels of the fast oscillators of the system, or the labels of the internal fields. 
To this end, we can now establish the main technical result of this paper. 

Theorem 11 Let G (A) and J^G (k) correspond to strictly Hurwitz stable open quantum 
systems (i.e., the A matrix of each system is strictly Hurwitz stable), and suppose that 
Sa + CiA~ 1 C*S- 1 — / and Sa — I are invertible. Then in the notation established above we 
have 

AFG (k) = J- AG (k) . 

The proof of the above theorem is given in the Appendix. Thus we establish that that if 
G (k) and J-"G (k) are strictly Hurwitz stable systems, and Sa + CiA~ l C* Si — I and Sa — I 
are invertible, the operation of adiabatic elimination of the oscillators in the network indeed 
commutes with the operation of taking the instantaneous feedback limit. For the systems in 



series example of section |IID| it can be seen that the strictly Hurwitz stable property holds 
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when A± and A2 are strictly Hurwitz stable, while for the beam splitter with an in-loop 



device of section II F it holds when en < 1. 



The requirement that G (k) and J-'G (k) be strictly Hurwitz is due to Remark 10. Note 
that the strict Hurwitz condition is not however necessary and that the limits may more 
generally commute whenever the kernel property of Y in Theorem [T] holds. 



V. CONCLUSION 



In this paper we have studied the question of commutativity of adiabatic elimination of 
oscillatory components and the operation of taking the instantaneous feedback limit in a 
quantum network with Markovian components. Provided some mild conditions are satis- 
fied, we answer the question in the affirmative by showing that adiabatic elimination can 
be viewed as a Schur complementation operation, thus putting it on the same footing as 
the instantaneous feedback limit, and subsequently proving the commutativity of successive 
Schur complementation. This result is important from a practical point of view because in 
practice it is much easier to obtain a simplified description of a quantum network by first 
obtaining simplified component models and then using them to obtain a description of the 
network rather than the converse operation of first forming the (possibly large) network and 
applying adiabatic elimination at the network level. Since we have shown that the order in 
which adiabatic elimination and the instantaneous feedback limit is taken is inconsequential, 
this justifies employing the former order of operations which is free of any concerns regard- 
ing the uniqueness of the resulting simplified network model in which the fast oscillatory 
components have been eliminated. 



VI. APPENDIX 

A. Proof of Theorem [T] 

Let us set 

M N = f) ® span I |n) : ^ rij = Nj , 

for N = 0, 1, 2, In particular, we have the direct sum of orthogonal subspaces f) <8> f) osc = 

©jv> Mjv. Let P s be orthogonal projection onto the "slow subspace" M = f) s = f) <8> C|0) osc 

29 



and let Pf — I — P s . Recall the hypothesis that ker(F) = f) s . We first have the following: 

Lemma 12 Under the hypothesis ker(F) = t) s , the subspaces M N are stable under Y N = 
Y\m n , and we have 

(PfYPf)- 1 = Om^Y- 1 . 

Moreover, let \5j) be the state where the j-th mode is in the first excited state and all others 
are in the vacuum, then (Yi) -1 ^\ <pj ® \Sj) = $ l ® I^j')- 

Proof. Stability and invertibility of Y N on M N follows directly from the specific form of 
Y N and the fact that Y has kernel space M Q . The relation (P f YP f ) 1 = ©^l*^ 1 follows 
from the direct sum decomposition. 

The remaining identity is easily checked from Y Y2j 4>j® \$j) = Ylji Aji<f>i<B)\5j) and setting 
this equal to ^ . <f>j <g> \5j) we deduce that 0/ = (A^ 1 )^ <f>j. m 

Corollary 13 ker(F*) = M„. 

Proof. By the preceding lemma we have that f)f = -P/f) ® f) OS c is stable under Y. Therefore, 
for any G M and ip G f) (g) t) osc we have that (0, Yip) = (0, YPfip) = 0. It follows that 
(Y*(p,ip) = (0, Yip) = for all ip G F) <g> f) OS c, thus Y*<p = for any G M and we conclude 
that M C ker(F*). We now need to show the converse that ker(F*) C M and we will do 
this by contradiction. To do this end, suppose that 3<p G Pft) <S> fw with <p ^ such that 
(Y*(p,ip) = for all ip G t) ® f) osc - It follows that (tp,Yip) = and therefore (tp,YP f ip) = 
for all -0 G f) <8> fW- But since f)f is stable under Y and Fj^ is invertible, it follows that 
if G t) s . But this contradicts the hypothesis that if is a non-zero element of ()/ and therefore 
we conclude that ker(F*) C M . This concludes the proof. ■ 

We now state a sufficient condition for ker(F) = M = ker(F*). Let us first recall the 
following definition 

Definition 14 A bounded Hilbert space operator A is strictly Hurwitz stable if 

Re{ip\Aip) < 0, for all ip ^ 0. 

Lemma 15 Let A jt G £>(f)) such that A = [A,-;] G £>(f) ® C m ) is strictly Hurwitz stable. 
Then the operator 

Y = ^A jl ®a* j a l (31) 
ji 

on f) <8> f) osc /ias kernel consisting of vectors of the form (g) |0) osc; where G f). 
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Proof. We see that for ip G \) eg) f) osc 

= (/ ® (A-J ® J ) ( J ® aO V>) = A * ® J ^> 

where ^ = (J ®bj)ip. We may decompose ^ = ^ (n) ® |n), where |n) is the orthonor- 
mal basis of number states for the oscillators and ipj (n) G t). Then 

n jl 

and, for each fixed n, we have Yljii^j ( n ) I Ajitpi (n)) < with equality if and only if the 
ipj (n) = since [Aji] is assumed to be strictly Hurwitz. In particular, if we assume that ip 
is in the kernel of Y then we deduce that ipj (n) = for each n and j — 1, • • • , m. It follows 
that ipj = (I Cg bj) ip — for each j = 1, • • • , m, and this implies that ip = </>® |0) osc for some 
</> G f) as required. ■ 

Note, however, that as we shall see below, for Theorem [T] to hold it is enough that 
ker(F) = M . 

Lemma 16 The operator S is unitary and K + K* + L*L = 0. 

Proof. We first show that I + CA^C* is invertible. Suppose that u G ker (J + CA~ 1 C*) 
u = -CA~ x C*u C*u = -C*CA- x C*u =>•(/ + C^CTT 1 ) C*w = 
(A + A~ 1 C*u = =► = C*w = 

so substituting C*u = into u = —C A~ l C*u we see that u — 0, therefore ker S 1 = 0. As 5 
is unitary, we have 

SS* = (I + CM^C) (/ + CA^C*) 

= I + CA' 1 {A + A* + A* -1 C* = I 



using the first of identities (12). Similarly S*S = I. 
Likewise we use ( 12 ) to show that 

K + K* + L*L = R- XA~ l Z + R* - Z*A*' X X* 

+ (G* - Z*A*- 1 C*) (G - CA- l Z) 

= - (X - G*C) A~ X Z - Z*A*~ X (X* - C*G - C*CA- 1 Z) 

= Z*A- l Z + Z*A*- 1 (Z + C*CA~ 1 Z) 

= Z*A*-\A + A* + C*C)A- 1 Z 

= 0. 
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Using the above results, we can now proceed to complete the proof of Theorem [TJ Let 
us first recall the main results from Bouten, van Handel, and SilberfarbP^. Let V (t, k) = 
U (t, k)*, then V satisfies the left QSDE (using a summation convention) 



dV{t, k) = V (t, k) {a (k) <g> dt + A (k) ® dBz(t) + 7, g> + foj - fy) ® dA jZ (*)} , 



where a (A;) = fc 2 a 2 + kaci + a = K (k)*, (3j (k) = kflij + /3 0j - = Lj (k)*, jj (k) = —Sf^Li 
and Eji = Sfj. Their results are stated for the left QSDE for the reason that it is easier to 
formulate the conditions for unbounded coefficients this way, however, the treatment is of 
course equivalent. 

We note that is then Y*, with kernel space M , and we denote its Moore-Penrose 
inverse by 0:2 (note that this Moore-Penrose inverse exists since Y* has the same form and 
properties as Y). The pre-limit coefficients satisfy Assumption 1 in the paper of Bouten, 
van Handel, and Silberfarb^ by construction. Assumption 2 of that work correspond to 
the identities a 2 a 2 = a 2 a 2 = Pf, <y 2 P s = 0, = and P s aiP s = 0: the last three 

are automatic since P s projects onto the ground state of the oscillator and in each case we 
encounter a 3 |0) osc = 0. The limit coefficients in Assumption 3 ofl^are then 



a = P s (ao - ai5 2 ai) P, = (R* - Z*A*~ l X*) ® |0)(0| osc = K* <g> |0)(0|, 
(3 = P s (f3 - a x &M P s = (G* - Z*A*- 1 C*) ® |0)(0| osc = L* ® |0}(0| os 
e = P s e (/ + $*a 2 ^) P s = S* (I + C*A*- l C*) ® |0)(0| osc = S* ® |0) (0| 
7 = -e/3* = -5*L®|0)(0| osc , 



with (S, L, K) as given in the statement of Theorem|lj These coefficients evidently satisfy the 
requirements of Assumption 3, namely to generate a unitary adapted Hudson-Parthasarathy 



equation on a common invariant domain in M , as was established in Lemma 16 
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B. Conditions for the Schur complements in Lemma [9] to be well-defined 
Lemma 17 If 



ker 



M BtB M B:C 
M c ,b M c , c 



C ker [M AJ3 M Afi \ 



lm 



M c ,a 



C im 



M b ,b M B , C 



ker M c ,c Q ker M B<C 

im Mc,b Q im Mc,c 
ker M B B C ker M C)B 
im Mb c C im M# b, 



(32) 

(33) 

(34) 
(35) 
(36) 
(37) 



tfiera t/ie Sc/mr complements (M/M c ,c)/(M/M c ,c)b,b and (M/M b ,b)/((M/M b ,b))c,c are 
all well-defined. 

Proof. Collecting the Schur complements used in the proof of Lemma [9] (successive 
complementation rule), we see that we have to show that 

M/M CtC , M/M BjB , (M/M c ,c)/(M/M c ,c)b,b 
Mauc,auc/M c ,c, Mauc,buc/M C) c, M Bu c,auc/M c ,c 



exist. To proceed, first note that, by Lemma 5l (32)-(37) imply that 



M/Mbuc,buc, M B uc, B uc/M BiB , Mb\jc,b\jc/Mc,c 

are well-defined. From ker Mb, b C ker Mc,b, we see that M B<B x = =>■ Mc )B x = 0. This 

M B B 



combined with condition (32) shows that M B b x — 



M, 



C,B 



x = M AB x = 0. 



Thus (39), given below, holds. Now, (33) implies that V x 3 y,z such that 



M BA 

Mq A 



X 



M BB y + M Bfi z 
M c , B y + M c , c z 



(38) 



But conditions (35) and (37) imply that 3 w,v such that Mc, B y = Mc,cv and M B ,cz 



M B B w. This together with (38) shows that imM B A Q imM B;B and imM Ct A Q imM^c 1 . 
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Combining this with (37) gives us (40), given below. By analogous arguments we can also 



establish (41) and (42) 



ker Mb b Q ker 



M A<B 
M c ,b 

im [ Mb a M B ,c] Q im Mb^b 

M A ,c 
M B ,c 

im[M c , A M c ,b] Q imM CjC 



ker Mcc Q ker 



(39) 
(40) 
(41) 
(42) 



From ([39]) to (J42j) it follows directly that M/M b ,b and M/M c ,c exist. Existence of 
Mauc,auc/Mc,c, M AuC ,buc/M c ,c and M B uc,Auc/Mc,c follows immediately from (Ell) and 



(42). 



Now we show that (M/Mc,c)/(M/Mc,c)b,b exists. We require 

ker(M B , B ~ M B ,cM^ c M c , B ) C ker(M A>B - M ac M^ c M c ,b) 
im(M B , A - M B ,cM- c M c , A ) C im(M B , B ~ M b ,cM^ c M ,b) 

Let v G im(M B>j4 - M B ,cM^ c M c , A ) be 

u = (M B)A - M B ,cMc <c M c , A )x 

M B ,ax 



[1 -M B aM ( 



c,c 



Mn A X 



for some vector x. Using (33) and noting that Mb,cM c c Mc,c 
Lemma [6]) leads to 



M B ,c (due to (34) and 



v = [l -M B nM. 



c,c 



M B , B y + M B , c z 
M c , B y + M C)C z 
= {M B , B - M b ,cM^ Mc,b)x 

which shows the required image space inclusion. 

To show that kei(M B ,B - Mb,cMc )C M c ,b) Q kev(M A>B - M ac M^ c M C) b) holds we 
choose some x G ker(M B B — Mb,c^c c^c,b) and see that 



M B , B x - M B ,cMc tC M a ,BX = [ M B , B M B , 



c 



X 



-Mc )C M c ,bx 
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which implies that 



x 



-M^ c M c , B x 



e ker[ M B)B M ByC ]■ However, M C)B x-M c ,cM cc M c , B x 



since Mc,cM cc Mc\ B = Mc, B (by (42) and Lemma [6j). It then follows that 



x 



G ker 



M BjB M B , C 
M C n M c ,c 



which by (32) implies 



x 



-XL. c M r B x 



G ker [Ma,b Ma,c) , an d therefore we deduce 



(Ma,b — Ma,cM c c Mc, B )x = which proves the required kernel space inclusion. 



We remark that the conditions (32)-(37) are not necessary, as is clear from Horn and 
Zhang 7 , page 42. 



C. Proof of Theorem El 



The model may initially be described by the set of coefficients g over the space 



ft = (f) © f)) © (C © £ c © Si) 



and we decompose this as 



ft = fi X © j} 2 © ^3 © f)4 

f) © (C © j?e) , Slow External 



where 

Si 

•f>2 = F) © ^ij Slow Internal 

j5 3 = F) © (C © £ c ) , Fast External 
#4 = F) © Fast Internal 

With respect to this decomposition, we decompose g into sub-blocks as 







(43) 



011 012 013 014 
021 022 023 024 
031 032 033 034 
,041 042 043 044 

The set of labels 3 = {1,2,3,4} can be split up into the slow labels {1,2} and the fast 
labels F = {3,4} = 5" as well as the external labels E = {1,3} and the internal labels 
/ = {2,4} = E'. 



Rx Afi X x 

G\ S[[ — / C[ 

Z x -C*Si A 
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To proceed, we must first establish that the generalized Schur complement is well-defined. 
Here we are ultimately retaining the "slow external" degrees of freedom (index 1) and 



eliminating the index sets {2, 3, 4}. To this end, We need to check that conditions (32)-(37) 



Let (x, y, z) T be an element of ker 



Then x, y satisfies (Sa — I)x + C{y 



are all satisfied. We begin with (]32|). 

Sa - I Q 
-C*S- X A 


and — C*S\X + Ay = 0, while z is arbitrary. Therefore, we have y = A~ 1 C*S- 1 x and 
(Sh + dA- x C*Si - I)x = 0. Since S a + CiA^CSi - I is invertible by hypothesis, we find 

s-i - 1 a 

that x = 0. It then follows that also y = 0, and we conclude that ker —C*S\ A 



consists of vectors of the form (0, 0, z) T . Clearly such vectors lie in the kernel of [ M\ X\ ] 

Sn-i a o 

and we conclude that ker — C*Si A ^ ker[ Mi Xi ]• 



Next, we check if for every given vector x there exist vectors y and z such that we have 
the equality 

n. nr q.. _ t r. 

(44) 







s-n - 1 a 




y 


Zix 




-C*Si A 




z 



In particular, this will be satisfied if the matrix 



Sii— I+CiA C*Si and A are invertible we see that this simply follows from the Banachiewicz 



Sa — I C\ 
~C*Si A 



is invertible. However, since 



formula. Therefore, for any vector x there indeed exist vectors y and z such that (44) holds 



invertible we also get 





G 1 




Sa 


a 







and we conclude that im 


Zi 


C im 


-C*Si 


A 





. Moreover, from the fact that A 





















ker 



im 



A 


-C*Si 




C ker[Ci 0], 



C im 



A 
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while from the invertibility of Sa — I we automatically have 



ker(5'ii — I) C ker 
im [ a ] C im(5ii - /) 



C*Si 




Therefore conditions (32)- (37) are satisfied, and the theorem now follows from Lemmata 17 
andU 
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